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We present cluster decay studies of 56 Ni* formed in heavy-ion collisions using different Skyrme 
forces. Our study reveals that different Skyrme forces do not alter the transfer structure of fractional 
yields significantly. The cluster decay half-lives of different clusters lies within ±10% for PCM and 
S ' ±15% for UFM. 
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I. INTRODUCTION 

Recently, a renewed interest has emerged in nuclear physics research. This includes low energy fusion process 
intermediate energy phenomena @ as well as cluster-decay and/or formation of super heavy nuclei [!, HJ . In the last 
one decade, several theoretical models have been employed in the literature to estimate the half-life times of various 
exotic cluster decays of radioactive nuclei. These outcome have also been compared with experimental data. Most 
of these models applied to study exotic cluster decay can be classified into two categories: In the first category, only 
barrier penetration probabilities are considered. Such models have been labeled as unified fission models (UFM) [5|- 
7-H ■ 0]. In the second category, clusters are assumed to be formed well before penetration. This is done by including the 
' preformation probability in the calculations. These models have been dubbed as preformed cluster models (PCM) 0- 
0\ , [lCj . In either of these approach, one needs complete knowledge of the potential. 

This problem is tackled in the literature in two different manners: One tries to adjust various parameters of model 
to known experimental data flll - [l4j . Alternatively, one starts from a basic fundamental approach free from such 
adjustable parameters [l5l - [2lj |. It remain to be seen how particular set of model parameters influence the cluster 



decay process. We plan to address this question in this paper. We shall work out the above problem with potential 



o 
o 



obtained from the Skyrme interactions. The Skyrme interactions are well used to describe the fusion process at low 
incident energies as well as in subthreshold, collective flow, and multifragmentation at intermediate energies. 

The Skyrme force is an effective interaction, which parameterizes the G-matrix by a zero ran ge, density and 
momentum dependent ansatz. The Skyrme force consists of two-body as well as three-body parts as |20| : 



^ i<j i<j<k 

Using a short-range expansion of the two-body interaction, the matrix elements in momentum space can be written 
as: 

(k | via I k') = t (1 + x P a ) + -ti (k 2 + k' 2 ) + t 2 k ■ k' + iW (oi + ct 2 ) • k x k', (2) 

where k and kl are the relative wave vectors of the nucleons. P G is spin exchange operator and a are Pauli spin 
matrices. To deal with such interaction, it is convenient to write the matrix elements in configuration space as: 

vxi = t (1 + x„P a ) S {f x - f 2 ) + [5 {fx - f 2 ) k 2 + k' 2 5 {fx - r 2 )] 

+t 2 k' ■8{fx-f 2 )k + iW {ax + 3 2 ) ■ k' x 8 {fx -f 2 )k, (3) 
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here k(= (Vi — V 2 )/2i) denotes the relative momentum operators acting on the right and fc'(= — (Vi — V2)/2t), 
acting on left, respectively. 

The three-body term of the Skyrme force can be written as: 



«i23 = *3<5 (n - r 2 ) 5 (f 2 - r 3 ) . 



(4) 



For the Hartree-Fock calculations of even-even nuclei, this force is shown to be equivalent to a two-body density 
dependent interaction: 



The above form, Eq. (JS|), provides a simple phenomenological representation of many body effects describing the way, 
in which the interaction between two nucleons is influenced by the presence of others. The Skyrme interaction is 
an approximate representation of the effective nucleon force which is valid only for the low relative momentum. In 
Eqs. @ to ([5]), we see several constants/parameters like t , ii, t 2l t 3 , xo, and Wq that need to be fitted. These 
parameters have been fitted by various authors from time to time to get better description of various ground state 
properties of nuclei [l7l [2p| - [22l | . A particular set comprising these parameters is known as Skyrme force. Till to-date, 
large number of Skyrme forces are available in the literature [22j . These different Skyrme forces constituting different 
equation of state at intermediate energies [231 l24j] . All the conventional (i.e. with the three body term replaced by a 
density dependent two body term), generalized (adjusting the effective mass m* and compressibility K) and modified 
Skyrme forces (adjusting the density parameter £3 to fit the spectra) are unified in a single form by Zhuo 125] as an 
extended Skyrme force. 



Our aim here is to study the role of various Skyrme forces and surface corrections in the exotic cluster decay process. 
This study is still missing in the literature. 

In recent years, there have been a number of experimental and theoretical studies [26M38I ] aimed at understanding 
the decay of light compound nucleus formed through heavy-ion reactions. In most of the reactions studied, whereas 
the general conclusion about the formation probability for the compound nucleus and characteristic features of its 
decay are debated in terms of either fusion-fission mechanism [26l [28l |3^ | , which may be considered as the emission 
of complex (or intermediate mass) fragments, or a deep inelastic (DI) orbiting po| mechanism behaviour. 

One of such system is the doubly magic 56 Ni, which is studied by using several entrance channels J_fO + 40 Ca, 
32 S + 24 Mg, 28 Si + 28 Si) and at different incident energies (1.5 to 2.2. times Coulomb barrier) [2^ - l36j . At these 
incident energies, the incident flux get trapped that results in the formation of compound nucleus, which is in addition 
to a significant large angle scattering cross-sections. For light masses (A<44), the compound nucleus decays by the 
emission of light particles and 7-rays. An experimental measure of this so called particle evaporation residue is the 
compound nucleus fusion cross-section. For heavier systems, such as 56 Ni, a significant decay strength to heavier 
fragments is also observed which could apparently not arise from a direct reaction mechanism because of large mass 
asymmetry differences between the entrance and exit channels. The measured angular distributions and energy spectra 
are consistent with fission like decays of the respective compound systems. 

The measured mass distribution for 56 Ni shows a preferential decays to channels comprising a-nuclei 16 0, 20 Ne, 



2 Mg and 28 Si, and their complimentary fragments [3J-l36|, independent of the entrance channel nuclei and centre- 
of-mass energy E crr ,. Such an a-structure is associated with the shell effects in the potential energy surface of the 
compound nucleus |35| . though these are almost zero at the compound nucleus excitation energies involved. Such an 
a-nucleus structure in the measured mass distribution of 56 Ni has its origin in the macroscopic energy [38[ . 

Cluster decay is studied for 56 Ni, when formed as an excited compound system in heavy-ion collisions. Since 56 Ni 
has negative Qout, and hence stable against both fission and cluster decay processes. However, if is is produced in 
heavy-ion reactions depending on the incident energy and angular momentum, the excited compound system could 
either fission, decay via cluster emissions or results in resonance phenomenon. The negative Q ou t is different for 
various exit channels and hence would decay only if it were produced with sufficient compound nucleus excitation 
energy E^, N (= E cm + Qi n ), to compensate for negative Q ou t, the deformation energy of the fragments Ed, their total 
kinetic energy (TKE) and the total excitation energy (TXE), in the exit channel as: 




(5) 




(6) 




E* CN =1 Qout I +E d + TKE + TXE. 



(7) 



(see Fig. 1, where Ed is neglected because the fragments are considered to be spherical). Here Qi n adds to the entrance 
channel kinetic energy E cm of the incoming nuclei in their ground states. 
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Section |TT] gives some details of the Skyrme energy density model and preformed cluster model and its simplification 
to unified fission model. Our calculations for the decay half-life times of 56 Ni* compound system and a discussion of 
the results are presented in Section IIIII Finally, the results are summarized in Section IIVI 



II. MODEL 



A. Skyrme Energy Density Model 

In the Skyrme Energy Density Model (SEDM) , the real part of interaction potential Vjy (R) is defined as difference 
between energy expectation value E of the whole system calculated at a finite distance R and at infinity [U [H[ . 

V N (r) = E (r) - E (oo) , (8) 

with 

' H(r)dr. CM 
In this formalism, the energy density functional H (r) read as; 

H(p,t,J) = A-T+iio[(l + ^o)p 2 -(xo + i)(p^+p^)] + |(ti+i 2 )pT 

+ g(*2 ~ tl)(p n T n + p p Tp) + -^(t 2 - iti)pV 2 p 
+ T^(3il + t 2 ){pn^ 2 Pn + /OpV 2 p p ) + ^t 3 p n ppp 

-iwoGoV • J + Pn V ■ J n + p P V • Jp). (10) 

Here p — p n + p p is the nucleon density taken to be two-parameter Fermi density and J = J n + J p is the spin density 
which was generalized by Puri et al. [17j . for spin- unsaturated nuclei. The remaining term is the kinetic energy density 
t = r n + t p . The Coulomb effects are neglected in the above energy density functional, but will be added explicitly. 
In Eq. fjlO[) . six parameters to, t\, t 2 , £3, xq, and Wo are fitted by different authors to obtain the best description 
of the various ground state properties for a large number of nuclei. As discussed in the introduction, these different 
parameterizations have been labeled as S, SI, SII, SIII etc.. 

The evaluation of kinetic energy density term was done within the Thomas-Fermi (TF) approximation which is a 
well known alternative to the Hartree-Fock method. As shown by various authors [4l[, the kinetic energy density r 
can be separated into volume term To and surface term plus reminder. In other words, 

T = T +TX + (11) 

In the first order approximation, one can limit to To term only. The volume term To in this approximation is given by 

3/3 5 



T0= 5\2> P ' 3 ' (12) 



The kinetic energy density r [4l|, |42[ , after including additional surface effects is 

2 



T = T + A^ (13) 

p 

here, A is a constant whose value has been a point of controversy and different authors have suggested different values, 
lying between 1/36 and 9/36. The above Thomas-Fermi approximation for t reduces the dependence of energy density 
H(r) to nucleon density p only. The exchange effects due to anti-symmetrization can be assimilated to reasonable 
extent when Eq. (1131) is used |17| . We apply the standard Fermi mass density distribution for nucleonic density: 

Pi{R) = ^ y, -oo<i?<oo (14) 



l + exp{^}' 
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The average central density poi given by [If 



POi 



3A, 



1 



(15) 



Roi and &i are, respectively, the half-density radii and surface diffuseness parameters taken from Refs.[lJ 
the details of the model, reader is referred to Ref. [l7j . 



For 



B. The Preformed Cluster Model 



For the cluster decay studies, we use the Preformed Cluster Model (PCM) This model, based on the quantum 

mechanical fragmentation theory [44T4471 ] . uses the decoupled approximation to 77- and i?-motions. The decay constant 
(A) in the PCM is defined as, 



K = v Q PP , 



or 7\ 



/2 



In 2 
~A~ 



(16) 



here Pq is the assault frequency with which the cluster hits the barrier, P is the probability of penetrating the barrier 
and Po is the preformation probability. Thus, in contrast to the unified fission models [5j-l2|, the two fragments in 
PCM are considered to be formed at a relative separation co-ordinate R before the penetration of the potential barrier 
with probability Po- The Schrodinger equation in terms of r\ and R coordinates as: 

H{r j ,R)i>(r ) ,R)=Eij(r,,R), (17) 

The above equation can be solved in a decoupled approximation [9| , for which the Hamiltonian takes the form: 

h 2 d 1 d h 2 d 1 d 



H = — 



^\fB~KR. dR VBrr. dR 



V( V )+V(R). 



(18) 



Since the potentials are calculated within the Strutinsky re-normalization procedure (V = VMacro + SU) by using an 
appropriate liquid drop model potential VMacro and asymmetric two center shell model for shell corrections SU, are 
nearly independent of the relative separation coordinate _R, R can be taken as a time independent parameter. For the 
Hamiltonian Eq. ljI5]l. the Schrodinger Eq. (fl~7| can be separated in two co-ordinates 77 and R as follows: 




and 



d 



1 



d 



2y/B^dRy/B^dR 



V(v) 



V(R) 



V>(i?) = E R ip(R), 



(19) 



(20) 



with ip(r), R) = ifj(r)) ip(R) and E = E v + E R . 

The fragmentation potential (or collective potential energy) V(rj), appearing in Eq. P^)l . is calculated as, 



rp2 \ 

V M acro{Ai,Zi) + SUi CXp ( -— 2 



Zi ■ Z 2 e 2 
R 



V N (R) + V t , 



(21) 



where the theoretical binding energies (V = VMacro + 5U) are taken from Moller et al. [48(]. The charges Zi in 
Eq. (|2ip arc fixed by minimizing the potential V(r/z), defined by Eq. (|2T|) without Vn{R) in r\z co-ordinates. The 
shell corrections 8U are considered to vanish exponentially for Eq N > 60 MeV, giving T = 1.5 MeV. At higher 
excitation energies, the shell corrections vanish completely and only the liquid drop part of energy is present. The 
additional attraction due to nuclear interaction potential Vn (R) is calculated within SEDM potential. The rotational 
energy due to angular momentum effects Ve (= h 2 £(£+l) /2/zi? 2 ) is not added here since its contribution to the structure 
yields is shown to be small for lighter systems [46]. The nuclear temperature T (in MeV), is related approximately to 
the excitation energy E^, N , as: 



E, 



CN 



-AT - T 



(in MeV). 



(22) 
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The kinetic energy part of the Hamiltonian in Eq. (fT9)) co mes through the mass parameter B nrt which is calculated 
using the classical mass parameter of Kroger and Scheid [49(, based on the hydrodynamical flow. The mass parameter 



B n „ reads as: 



with 



R, 



2R„ 



B 



AmR min 



vv 



1 



1 ^ c 
cos »i v Ri 



1 



R c 

COS#2 V 



(23) 



(24) 



2R„ 



[(1 - COS0i)(i?i - Rc) + (1 - COS0 2 )(i? 2 - R c )} , 



(25) 



7T_R C Rrnin j 



i? r = 0.4i? 2 



(26) 



and i>t = ui + «2, is the total conserved volume. 

Solving Eq. (19) numerically, | ip(rj) \ 2 gives the probability of finding the mass fragmentation r\ at a fixed position 
R, on the decay path. Normalizing and scaling | ip(rj) | 2 to give the fractional mass yield for each fragment in the 
ground state decay as: 



The nuclear temperature effects in Eq. (I27p are also included through a Boltzmann-like function, 



i 2 = X! i i 2 exp 



T 



(27) 



(28) 



For R- motion, instead of solving the stationary Schrodinger Eq. (|20p . the WKB action integral was solved for the 
penetration probability P [5(| ■ For each rj- value, the potential V(R) is calculated by using SEDM for R > i?^, with 
Rd = Rmin + Ai? and for i? < i?^, it is parameterized simply as a polynomial of degree two in R: 



V(R) 



I Qout I +ai(R - i? ) + a 2 (R - R Q ) 2 for R Q <R< R d , 
V N (R) + Z x - Z 2 e 2 /R for R > R d , 



(29) 



where R is the parent nucleus radius and AR is chosen for smooth matching between the real potential and the 
parameterized potential (with second-order polynomial in R). A typical scattering potential, calculated by using 
Eq. (29) is shown in Fig. 1, with tunneling paths and the characteristic quantities also marked. Here we choose the 
first (inner) turning point R a at the minimum configuration i.e. R a — Rmin (corresponding to V m i n ) with potential 
at this i? a -value as V(R a = Rmin) = V m in (displayed in Fig. 1) and the outer turning point Rb to give the Q e //-value 
of the reaction (Q e ff =| Qout | +TKE) i.e. V(Rb) = Q e ff- This means that the transmission probability P with the 
de-excitation probability, Wi = exp(— WSj) taken as unity, can be written as: 



P = P l Pb, 

where P, and Pj, are calculated by using WKB approximation, as: 

Ri 

Pi = exp 



{2fx[V(R) - V(R l )]} 1 ' 2 dR 



(30) 



(31) 



and 



exp 



{2n[V(R)-Q eff ]} 1/2 dR 



(32) 



6 



here R a and Rb are, respectively, the first and second turning points. This means that the tunneling begins at 
R = Ra (= Rmin) an d terminates at R = Rb, with V(Rb) — Qeff- The integrals of Eqs. (31) and (32) are solved 
analytically by parameterizing the above calculated potential V(R). 

The assault frequency or the barrier impinging frequency vq in Eq. (16), is given simply as, 

where E% = ^Qeff is the kinetic energy of the emitted cluster, with Q e ff shared between the two fragments and 

/j, = m{ a ) is the reduced mass. 

The PCM can be simplified to unified fission model (UFM), if preformation probability Pq = 1 and the penetration 
path is straight to Q e //-value. 



III. RESULTS AND DISCUSSIONS 



The calculations are made in two steps: In the first steps, we studied the role of different Skyrme forces in the 
cluster decay of 56 Ni* and in the second step, effect of surface correction term A is analyzed. 

Fig. 1 shows the characteristic scattering potential for the cluster decay of 56 Ni* into 16 + 40 Ca channel as an 
illustrative example. In the exit channel for the compound nucleus to decay, the compound nucleus excitation energy 
Eq N goes in compensating the negative Q ou t, the total excitation energy TXE and total kinetic energy TKE of the 
two outgoing fragments as the effective Q- value (i.e. TKE = Qeff in the cluster decay process). In addition, we plot 
the penetration paths for PCM and UFM. For PCM, we begin the penetration path at R a — R m in with potential 
at this i? a -value as V(R a = Rmin) = V m in and ends at R = Rb, corresponding to V(R = Rb) = Qeff, whereas 
for UFM, we begin at R a and end at Rb both corresponding to V(R a ) = V{Rb) = Qeff- We have chosen only the 
case of different Q e ff (listed in Table 1), for different cluster decay products to satisfy the arbitrarily chosen relation 
Q eff = 0.4(28- | Q out |) MeV, as it is more realistic p| . 



A. Role of Different Skyrme Forces 

Figs. 2(a) and (b) shows the fragmentation potential V(r]) and fractional yield at R = R m i n with V(R m i n ) = V m i n . 
The classical hydrodynamical mass parameter B vv of Kroger and Scheid |49j used in the calculation of preformation 
probability. The fractional yields are calculated within PCM at T = 3.0 MeV using different Skyrme forces for 56 Ni*. 
From the figure, we observe that different Skyrme forces do not alter the transfer structure of fractional yields. The 
Skyrme force parameters have marginal role to play. Some variations in the absolute values are however visible [5ll | . 
The fine structure is not at all disturbed for different sets of Skyrme forces. 

The results for the cluster decay half-lives in 56 Ni* are quantified by the following quantity as: 

n rr lo, (iogri /2 r-(iog?\/2) s/// im .... 

[log T 1/2 ] % = ( logri/2) sm x 100 ' ( 34 ) 

where i stands for different sets of Skyrme force parameters and SIII for one set of Skyrme force parameters, which 
is widely used. Here, the strength parameter of surface correction is taken as zero (i.e. A = 0). 

In Fig. 3(a) and (b), we display the quantified results using Eq. ([34|) for logT^ within PCM and UFM models as 
a function of cluster mass A 2 . The role of temperature T (or excitation energy Eq N ) enters only in the PCM via 
preformation probability Pq. These variation in the cluster decay half-lives for different clusters lies within ±10% for 
PCM and ±15% for UFM. This amount is significant once we understand cluster decay probabilities can be measured 
with great accuracy in the literature. 



B. Role of Strength Parameter of Surface Correction (A) 



The effect of different A- values for the heavy- ion nuclear potential is analyzed in Refs. [17|, l52j, suggesting that 
different A- value, can alter the depth of the nuclear potential Vjv significantly. In Ref. jl7[, it was shown that the 
barrier heights gets lowered whereas the fusion barrier position shifts outward where stronger role of A is taken into 
account. The effect of this strength parameter A for additional surface effects in the decay calculations has yet not been 
studied in the literature. In this subsection, we plan to study the effect of strength parameter of surface correction 
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on cluster decay half- lives by taking different A- values (equal to 0, 1/36, 2/36, 3/36, 4/36, and 5/36) in SEDM for 
the compound system 56 Ni*. 

In Fig. 4, the scattering potential for different values of surface correction factor A is plotted as a function of 
internuclear distance R. One observes from the figure that variation in the A-value changes the interior part of the 
scattering potential thereby changing the penetration probability. 

In Fig 5(a) and (b), we show the fragmentation potential V(?7)and fractional mass distribution yield at R = R m i n 
with V(R m in) — V m in- The fractional yields are calculated within PCM at T = 3.0 MeV using different values of 
surface correction factor for 56 Ni*. From figure, we observe that different values of A changes the fractional yield to 
large extent but do not alter its transfer structure. The fine structure is not at all disturbed for different values of 
surface correction factor. 

The results for the cluster decay half-lives in 56 Ni* are quantified by the following quantity as: 

Pi t 1 of (logT 1/2 )^(logT 1/2 )^ 

[logT 1/2 ] % = (i ogri/2 )A=o X 100 < ( 35 ) 

where i stands for different A-values of the strength parameter of surface correction. Skyrmc force SIII is employed for 
these calculations. In Fig. 6, we display the quantified results using Eq. (35) for the percentage variation of logTi/ 2 
within PCM and UFM as a function of cluster mass A 2 . The variation in the cluster decay half- lives for different 
clusters lies within ±10% for both PCM and UFM. Together with the effect of different Skyrme forces, one can see 
that the net effect of different Skyrme forces as well as surface corrections has sizable effect on the cluster decay 
half-life times. 



IV. SUMMARY 



We here reported the role of different Skyrme forces as well as surface corrections in the cluster decay constant 
calculations. Our studies revealed that the effect of different Skyrme forces on the cluster decay half-life times is 
about ±15%, whereas it is ±10% in the case of surface corrections. 

This work was supported by a research grant from the Department of Atomic Energy, Government of India. 
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TABLE I: The calculated characteristic quantities for cluster decay of 56 Ni* compound system for fragment masses A2 > 16, 
with excitation energies E* — Q e ff+ \ Qout |. 



Cluster + 


| Qout 1 


Qeff 


E* 


Daughter 


(MeV) 


(MeV) (MeV) 


ib O + 4U Ca 


-14.12 


5.55 


19.67 


^Ne + a8 Ar 


-22.23 


2.31 


24.54 


^ u Ne + ' w Ar 


-17.12 


4.35 


21.47 


^Mg + a4 S 


-24.58 


1.37 


25.95 


M Mg + 


-16.57 


4.57 


21.14 


*Si + 3,) Si 


-23.57 


1.77 


25.34 


' M Si + ' M Si 


-12.20 


6.32 


18.52 



10 




^5r 10 15 20 25 30 35 40 45 50 
** a R (fm) 



FIG. 1: The scattering potential V(-R) (in MeV) for cluster decay of 56 Ni* into 16 + 40 Ca channel using Skyrme force SIII, 
with A = 0. The distribution of compound nucleus excitation energy Eqn a t both the initial (R = Ro) and asymptotic (R — i> oo) 
stages and Q-values are shown. The decay path for both PCM and UFM models is also displayed. 
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FIG. 2: (a) The fragmentation potential V(rj) and (b) calculated fission mass distribution yield with different Skyrme forces 
at T = 3.0 MeV. 
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FIG. 3: Percentage variation of logT^ for different Skyrme forces w.r.t. SIII force. 
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FIG. 4: Same as Fig 1, but for different values of surface correction factor (A). 
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FIG. 5: Same as Fig 2, but for different values of surface correction factor (A). 



13 



56 Ni* 



A 

O 



(a) 



o 

A 



A 



PCM 



O 
A 
V 



o 

A 



V 

O 



O 
< 



Surface Correction (X) 

O A=l/36 A X=2/36 
V X=3/36 O X=4/36 
< X=5/36 



o 

A 



O 
A 



V 

O 



(b) 
_i 



o 

A 



15 



10 



^ 
o 

-5 " 

-10 

-15 



16 18 20 22 24 26 28 16 18 20 22 24 26 28 

Fragment Mass A 2 



FIG. 6: Same as Fig 3, but for different values of strength parameter of surface correction (A =fc 0) w.r.t. surface correction 
(A = 0). 



